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The RICIS Concept 


The University of Houston-Clear Lake established the Research Institute for 
Computing and Information Systems (RICIS) in 1986 to encourage the NASA 
Johnson Space Center (JSC) and local industry to actively support research 
in the computing and information sciences. As part of this endeavor, UHCL 
proposed a partnership with JSC to jointly define and manage an integrated 
program of research in advanced data processing technology needed for JSC’s 
main missions, including administrative, engineering and science responsi- 
bilities. JSC agreed and entered into a continuing cooperative agreement 
with UHCL beginning in May 1986, to jointly plan and execute such research 
through RICIS. Additionally, under Cooperative Agreement NCC 9-16, 
computing and educational facilities are shared by the two institutions to 
conduct the research. 

The UHCL/RICIS mission is to conduct, coordinate, and disseminate research 
and professional level education in computing and information systems to 
serve the needs of the government, industry, community and academia. 
RICIS combines resources of UHCL and its gateway affiliates to research and 
develop materials, prototypes and publications on topics of mutual interest 
to its sponsors and researchers. Within UHCL, the mission is being 
implemented through interdisciplinary involvement of faculty and students 
from each of the four schools: Business and Public Administration, Educa- 
tion, Human Sciences and Humanities, and Natural and Applied Sciences. 
RICIS also collaborates with industry in a companion program. This program 
is focused on serving the research and advanced development needs of 
industry. 

Moreover, UHCL established relationships with other universities and re- 
search organizations, having common research interests, to provide addi- 
tional sources of expertise to conduct needed research. For example, UHCL 
has entered into a special partnership with Texas A&M University to help 
oversee RICIS research an 4 education programs, while other research 
organizations are involved via the “gateway* concept 

A major role of RICIS then is to find the best match of sponsors, researchers 
and research objectives to advance knowledge in the computing and informa- 
tion sciences. RICIS, working jointly with its sponsors, advises on research 
needs, recommends principals for conducting the research, provides tech- 
nical and administrative support to coordinate the research and integrates 
technical results into the goals of UHCL, NASA/JSC and industry. 
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1.0 Introduction : 


As part of the RICIS activity, the reinforcement learning techniques developed at 
Ames Research Center are being applied to proximity and docking operations 
using the Shuttle and Solar Max satellite simulation. This activity is carried out 
in the software technology laboratory utilizing the Orbital Operations Simulator 
(OOS). 

This interim report provides the status of the project after two months of activities 
and outlines the future plans. Our technical progress is summarized in section 2 
and future plans are described in section 3. 


2.0 Summary of Technical Progress : 

Considerable technical progress, as described below, has been achieved 
during these two months of activities in the project. 

1. The source code received from Ames Research Center contained the attitude 
controller however, the equations of motion still pertained to the inverted 
pendulum motion. In order to reflect the shuttle motion in the orbital 
environment, the source code was modified to interface with appropriate 
routines in the OOS. First, the onorbit digital auto pilot was modified to call the 
learn cycle in its normal calling sequence. Second, a special OOS structure 
was created to hold all relevant parameters for input and output. Third, an 
initialization routine was created to initialize the parameters of the learning 
algorithms. And, finally, the output of the learn cycle was interpreted and 
converted to rotational hand controller command. The source code is provided 
in Appendix A. 

2. An attitude hold test case was set up with the shuttle in its normal orbit during 
its mission. The pitch attitude of the shuttle was set at 45 degrees in the local 
vertical local horizontal coordinate system. Because of the gravity gradient 
torques, the shuttle pitch rate starts increasing resulting in an increase in the 
pitch attitude. As soon as the attitude error (difference between the desired 
attitude and the current attitude ) goes beyond the deadband limit, the fuzzy 
controller will initiate jet firing. Two neural networks are learning during each 
cycle and adjusts the parameters as required. The test case was setup to run up 
to 1000 seconds so there will be at least three attitude firings. For debug and 
analysis purposes, we initially performed the test case for only 30 seconds. 
Later, we extended the test case to 1000 and 100,000 seconds to understand 
the learning process and evaluate the performance of the fuzzy controller and 
neural networks. Plots for the 1000 seconds test case are included in Appendix 
B. This is one of the typical tests performed. So far we have performed over 10 
test cases with modifications in the source code and changes in the learning 
parameters. Results and our findings are described separately in this section. 



3. During the debug phase, we found several implementation errors. We had to 
change our structure, interface code and membership functions definitions. Our 
initial test results showed that the algorithm is holding the rate but not the 
attitude. We checked for all possible implementation errors, and when we were 
satisfied that there are no critical implementation errors, we did the test case 
again. 

4. When we found that the results are the same (i.e. no attitude hold), we started 
analyzing the Tsukamoto's defuzzification method used in the algorithm. In this 
method, only one side of the triangle is used in computing the output action. In 
our attitude controller, we have used the max-dot and max-min methods to 
compute the output, and thus we compared the two methods for any differences. 
We found that the Tsukamoto's method computes the output with only half the 
value of the max-dot method. Thus, we changed the interpretation of the 
learning output for generating the rotational hand controller command. Results 
of our test indicated that the controller now holds attitude as well as rate. The 
analysis was involved and we learned a lot about the algorithm. 

5. Since our test was successful in holding the rate and attitude, we started 
analyzing the learning behavior neural networks. Our results indicate that the 
parameter T that relates to the firing strength varies in the same manner for all 
rules. The 'd* parameters also show the same variations for all rules. We have 
discussed the problem with Dr. Berenji at Ames and sent him the source listing 
as well as plots so that he can compare it with his results and tell us if we are 
doing something wrong. 

6. We have analyzed the effect of the bias term x(2) and have concluded that 
the its most appropriate value is 0.5. We performed tests with bias values of 0.0 
and 1 .0 and results were really bad. For these values the controller does not 
perform its functions and neural networks do not learn. 

7. We were invited by the SPIE fuzzy logic and neural network committee to 
present our initial results at the Orlando conference and special workshop on 
Fuzzy-Neuro control. We presented our preliminary results at this workshop on 
April 23, 1992. 


3.0 Future Plans : 

We plan to continue testing of the fuzzy learning algorithms utilizing the attitude 
hold test case in the orbital operations simulator. Emphasis during this activity is 
to ensure that the algorithms perform properly and learning by the neural 
networks in ARIC architecture is achieved in a satisfactory manner. Then we will 
switch to other attitude control tests such as attitude maneuver, rate hold and 
rate maneuver. We will perform these tests with proper perturbations typically 
present during the orbital operations. 

As soon as we complete testing the attitude controller with learning, we will 
implement the translational controller in the simulation and perform proximity 



operations test cases. Currently we plan to perform v-bar, r-bar, fly around and 
station-keeping test cases as we have performed these test cases to check out 
designs of our translational controller. Finally, we will set up a test case that will 
simulate docking operations. In this test case, the shuttle will approach the solar 
max satellite from 50 feet to 2 feet and hold the relative orientation for a 
specified time at the final distance so that the grappling task can be performed. 


4. Issues and Concerns : 

At this time the project is on schedule and there are no issues or concerns to 
report in this interim report. 



Appendix A. Source Listing of Fuzzy Learning Modules 
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#include <stdio.h> 

#include <math.h> 

#include <sys/types.h> 

r EXTERNAL DATA STRUCTURE DEFINITION 7 
#include " . Jo rb_f u zzy/l e ar n_cy cl e . h " 
int counter=0; 

#define max(x,y) ( (x >= y) ? x : y ) 

#define min(x,y) ( (x < y) ? x : y ) 

#define Gamma 0.9 
#define Beta 0.2 
#define Beta_h 0.05 

#define Rho 1 .0 

#define Rho_h 0.2 

#define Rhol 2.0 

#define Rho_h1 0.4 
extern double sgn(); 
extern double exp(); 
extern double rnd(); 

learn_cycle(L) 

LEARN_CYCLE * L; r IN : 7 

{ 

int i.il, j, k; 

double match () f calculate_z_array(); 

double temp ; 


r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
r 11 March 1 992 - Alter scaling 7 

r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
r L->x[0] = L->Phi/20.0; 7 

r L->x[1 ] = L->Phi_dot/4.0; 7 

r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
r 08 April 1 992 - Alter scaling 7 

r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
r L->x[0] = L->Phi/1 0.0; 7 

r L->x[1 ] = L->Phi_dot/2.0; 7 

r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
L->x[0] = L->Phi ; 

L->x[1] = L->Phi_dot * 10.0 ; 

r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
r 08 April 1992 - Alter Bias 7 

r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
r L->x[2] = 1.00; 7 
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r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 

r 15 April 1992- Alter Bias 7 

r CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 7 
L->x[2] = 0.50 ; 


L->failure = 0; 

r Set up and evaluate the failure criteria 7 
if ( (fabs(L->Phi) > 0.7) || (fabs(L->Phi_dot) > 0.07) ) 
L->failure = -1 .; 

r output: state evaluation 7 

for (i = 0; i < 31; i++) 

{ 

L->sum = 0.0; 
for(j = 0; j < 3; j++) 

L->sum += L->a[i*3+j] * L->x_old[j]; 

L->y[i] = 1 .0 / (1 .0 + exp(-1 .0 * L->sum)); 

} 

L->sum = 0.0; 
for(i = 0; i < 31 ; i++) 

L->sum += L->c[i] * L->y[i] ; 

} 

L->sum1 = 0.0; 
for ( j = 0;j < 3; ]++) 

^ L->sum1 += L->b[j] * L->x_old[j]; 

L->v = L->sum + L->sum1 ; 

r output: action 7 
for(i = 0; i < 31 ; i++) 

{ 

ii=i; 

L->w[i] = match (il ,L); 

L->zl [i] = calculate_z_array(i1 ,L); 

} 

L->sum1 = 0.0; 

L->denom = 0.0; 
for(i = 0; i < 31 ; i++) 
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{ 

L->sum1 += L->w[i] * L->z1[i] *L->f[i] ; 
L->denom += L->w[i]*L->f[i] ; 


} 

L->push = L->sum1 / L->denom; 

r output: action computations completed 7 
for(i = 0; i < 31 ; i++) 

{ 


L->sum = 0.0; 

for (j = 0; j < 3; j++) 

L->sum += L->d[i*3+j] * L->x_oldQ]; 

L->z[i] = 1 .0 / (1 .0 + exp(-1 .0 * L->sum)); 

L->sum2 = 0.0; 

L->sum3 = 0.0; 

for(i = 0; i < 3; i++) 

L->sum2 += L->e[i] * L->x_old[i]; 

for (i=0;i < 31 ; i++) 

L->sum3 += L->f[i] * L->z[i]; 

r ccccccccccccccccccccccccccccccccccccccccccccccccccc 

C 7 

r 4 May 1992 - Normalize sum3 7 

r L->sum4 = L->sum3 +L->sum2; 7 

rccccccccccccccccccccccccccccccccccccccccccccccccccc 

C 7 

L->sum4 = ( L->sum3 / 31 .0 ) + L->sum2; 

L->p = 1 .0 / (1 .0 + exp(-1 .0 * L->sum4)); 

rccccccccccccccccccccccccccccccccccccccccccccccccccc 

C 7 

r 1 5 April 1 992 - Use temp variable - not push 7 

rccccccccccccccccccccccccccccccccccccccccccccccccccc 

C 7 

r L->push = (rnd() <= L->p) ? L->push : -L->push; 7 
r L->unusualness = (L->push > 0) ? 1 .0 - L->p : -L->p; 7 
temp = (rnd() <= L->p) ? L->push : -L->push; 

L->unusualness = (temp > 0) ? 1 .0 - L->p : -L->p; 
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r using new input values and unmodified weights. 7 
I* Use y_new and v_new so not to destroy y and v. 7 
for(i = 0; i < 31 ; i++) 

{ 

L->sum = 0.0; 
for(j = 0; j < 3; j++) 

L->sum += L->a[i*3+j] * L->x[fl; 

L->y_new[i] - 1 .0 / (1 .0 + exp(-1 .0 * L->sum)); 

} 

L->sum « 0.0; 

L->sum1 = 0.0; 

L->sum2 = 0.0; 

for(j = 0; j < 3; j++) 

L->sum1 += L->b[j] * L->x[fl; 

for(i = 0; i < 31 ; i++) 

L->sum2 += L->c[i] * L->y_new[i]; 

L->sum = L->sum1 + L->sum2; 

L->v_new = L->sum; 

r action evaluation 7 
if (L->failure) 

L->r_hat = L->failure - L->v; 
else 

L->r_hat * L->failure + Gamma * L->v_new - L->v; 
r modification and update to parameters 7 
for(i = 0; i < 31 ; i++) 

L->factor1 = Beta_h * L->r_hat * L->y[i] * (1.0 - L->y[i]) * sgn(L->c[i]); 
L->c[i] += Beta * L->r_hat * L->y[i]; 

for(j - 0; j < 3; ]++) 

L->a[i*3+j] += L->factor1 * L->x_old[j]; 

} 

} 


for(i = 0; i < 3; i++) 
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L->b[i] += Beta * L->r_hat * L->x_old[i]; 

for(i = 0; i < 31 ; i++) 

{ 

L->factor2 = Rho_h * L->r_hat * L->z[i] * (1.0 - L->z[i]) * sgn(L->f[i])*L->unusualnes 

s; 

for(j = 0; j < 3; j++) 

L->d[i*3+j] += L->factor2 * L->x_oldO]; 

} 


for(i = 0; i < 31 ; i++) 

{ 

L->f[i] += Rho * L->r_hat * L->unusualness * L->z[i]; 
for(i = 0; i < 3; i++) 

{ 

L->e[i] += Rho * L->r_hat * L->unusualness * L->x_old[i]; 


L->x_old[0] = L->x[0]; 
L->x_old[lj « L->x[lj; 
L->x_old[2] = L->x[21; 



double sgn(x) 
double x; 

{ 

if (x < 0.0) 
return (-1.0); 
else if (x > 0.0) 
return (1.0); 
else 

return (0.0); 

} 


!* zero_one function returns 0 for negative numbers 

1 for values > 1 

x for values between 0 and 1 V 
double zero_one(x) 
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double x; 

{ 

if (x < 0) return (0.0); 
else if (x > 1) return (1 .0); 
else return (x); 

} 


r********** Membership Function for Phi ******************* 

double nbl(x) 
double x; 

return(min( max((-2-x)/3 , 0.0 ), 1.0 )); 

} 

double nml(x) 
double x; 

^ if (x <= -2.5) return (min( max(( x+4)/1 .5 , 0.0 ), 1 .0 )); 
else return (min( max(( -x-1)/1.5 , 0.0 ), 1.0 )); 

} 

double nsl (x) 
double x; 

if (x <= -1 .0) return (min( max( x+2.0 , 0.0 ), 1 .0 )); 
else return (min( max( -x , 0.0 ), 1.0 )); 

} 

double zol (x) 
double x; 

if (x <= 0) return (min( max( x+1 , 0.0 ), 1 .0 )); 
else return (min( max( -x+1 , 0.0 ), 1.0 )); 

} 

double psl (x) 
double x; 

if (x <= 1 ) return (min( max( x , 0.0 ), 1 .0 )); 
else return (min( max( -x+2 , 0.0 ), 1.0 )); 

} 

double pml(x) 
double x; 

^ if (x <= 2.5) return (min( max(( x-1 )/1 .5 , 0.0 ), 1 .0 )); 
else return (min( max(( -x+4)/1 .5 , 0.0 ), 1 .0 )); 
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} 

double pbl (x) 
double x; 

{ 

return(min( max(( x-2)/3 , 0.0 ), 1.0 )); 

} 

/********•********** Membership Functions ********* 

double nb2(x) 
double x; 

{ 

return(min( max((-.2-x)/.3 , 0.0 ), 1.0 )); 

} 

double nm2(x) 
double x; 

{ 

if (x <= -.25) return (min( max(( x+.4)/.15 , 0.0 ), 1 .0 )); 
else return (min( max(( -x-.l )/.1 5 , 0.0 ), 1 .0 )); 

} 

double ns2(x) 
double x; 

{ 

if (x <= -.1) return (min( max(( x+.2)/.1 , 0.0 ), 1 .0 )); 
else return (min( max(( -x/.l) , 0.0 ), 1.0 )); 

} 

double zo2(x) 
double x; 

{ 

if (x <= 0) return (min( max(( x+.l )/.1 , 0.0 ), 1 .0 )); 
else return (min( max(( -x+.1)/.1 , 0.0 ), 1.0 )); 

} 

double ps2(x) 
double x; 

{ 

if (x <= .1 ) return (min( max(( x/.l ) , 0.0 ), 1.0)); 
else return (min( max(( -x+.2)/.1 , 0.0 ), 1.0 )); 

} 

double pm2(x) 
double x; 

{ 

if (x <= .25) return (min( max(( x-.l )/.1 5 , 0.0 ), 1 .0 )); 
else return (min( max(( -x+.4)/.15 , 0.0 ), 1.0 )); 

} 


double pb2(x) 
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double x; 

return(min( max(( x-.2)/.3 , 0.0 ), 1.0 )); 

} 

^***************** Defu zz ification Process with Accel Membership Functions 

double nm3(x) 
double x; 

{ 

return(-2-x); 

} 

double ns3{x) 
double x; 

{ 

return(-2*x); 

} 

double zo3(x) 
double x; 

{ 

return(O.O); 

} 

double ps3(x) 
double x; 

{ 

return(2*x); 

} 

double pm3(x) 
double x; 

{ 

return(2+x); 

} 

double match(i.L) 
int i; 

LEARN_CYCLE *L; TIN : 7 

{ 

double temp; 
switch (i) { 


case 0: 


"* + 
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temp =min( zero_one(nb1(L->Phi)*L->d[i*3+0]), zero_one(zo2(L->Phi_dot)*L->d[i*3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp= min(temp,zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
case 1: 

temp=min( zero_one(nb1 (L->Phi)*L->d[i*3+0]) t zero_one(ns2(L->Phi_dot)*L->d[i*3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp,zero_one(L->x[2]*L->d[i*3+2])); 

V 

return(temp); 
case 2: 

temp=min( zero_one(nb1 (L->Phi)*L->d[i*3+0]) t zero_one(nm2(L->Phi_dot)*L->d[i*3+ 

1])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
case 3* 

temp=min( zero_one(nb1 (L->Phi)*L->d[i*3+0]), zero_one(nb2(L->Phi_dot)*L->d[i*3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
case 4: 

temp=min( zero_one(nm1 (L->Phi)*L->d[i*3+0]), zero_one(zo2(L->Phi_dot)*L->d[i*3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])) ; 

7 

return(temp); 
case 5: 

temp=min( zero_one(nm1 (L->Phi)*L->d[i*3+0]) t zero_one(ns2(L->Phi_dot)*L->d[i*3+ 

1 ])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
case 6: 

temp=min( zero_one(nm1 (L->Phi)*L->d[i*3+0J), zero_one(nm2(L->Phi_dot)*L->d[i*3 

i])); 

REMOVE "BIAS" REFERENCE 
temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 
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7 

return(temp); 

“temp-minf zero_one(nm1 (L->Phi)*L->dp*3 + 0]). zero_one(nb2(L->Phi 

-my 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 

Ca tlmp-min( zero_one(ns1 (L->Phi)*L->dp-3 + 0J). zero_one(zo2(L->Phi. 
11 ))' 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2] L->d[i 3+2])), 

7 

return(temp); 

“Imp-minf zero_one(ns1 (L->Phi)*L->d[i*3 + 0]). zero_one(ns2(L->Phi. 
11 ))' 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 

C ^mp-min( zero_one(ns1 (L^PhiTU.dO-3+0]). zero_one(nm2(L->Ph 
11 \\. 

T REMOVE "BIAS" REFERENCE Hn *o. oi\v 

temp=min(temp, zero_one(L->x[2] L->dp 3+2])), 

7 

return(temp); 

Templmint zero_one(ns1 (L->Phi)*L->dp*3 + 0]). zero_one(nb2(L->Ph 

1 ])); 


r REMOVE "BIAS” REFERENCE 

temD=min(temp, zero_one(L->x[2] L->d[i 3+2])), 


7 


temp=min(temp 
return(temp); 

C3Se ’j-mint zero_one(zo1(L->Phi)*L->dp*3+0]), zero_one(pb2(L->Pf 


temp= 


ID): 


r REMOVE "BIAS" REFERENCE 

temD=min(temp, zero_one(L->x[2] L->d[i 3+2])), 


7 


temp=min(temp 
return(temp); 

Temp-minf zero_one(zol(L->Phi)*L->dp*3 + 0]), zero_one(pm2(L->P 


dot)*L->d[i*3+ 


dot)*L->d[i*3+ 


dot)*L->d[i*3+ 


i_dot)*L->d[i*3+ 


i_dot)*L->d[i*3+ 


i_dot)*L->d[i*3+ 


i_dot)*L->d[i*3+ 
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i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 

temp=min( zero_one(zo1 (L->Phi)*L->d[i*3+0]), zero_one(ps2(L->Phi_dot) L->d[i 3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
case 15* 

temp=min( zero_one(zo1 (L->Phi)*L->d[i*3+0]). zero_one(zo2(L->Phi_dot)*L->d[i-3 + 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 

case ^0* 

temp=min( zero_one(zo1 (L->Phi)*L->d[r3+0]) f zero_one(ns2(L->Phi_dot)*L->d[i*3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 

C ^emp-min( zero_one(zo1 (L->Phi)*L->d[i*3+0]), zero_one(nm2(L->Phi_dot)*L->d[i*3+ 

1])): 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i 3+2])); 

7 

return(temp); 

C ^mp=min( zero_one(zo1 (L->Phi)*L->d[i*3+0]), zero_one(nb2(L->Phi_dot)*L->d[i*3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 

“temp=min( zero_one(ps1 (L->Phi)'L->d[i*3+0]), zero_one{pb2(L->Phi_dol)*L->d[i*3 + 
11 ))* 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 


7 
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return(temp); 


^4° mint zero_one(ps1 (L->Phi)"L->dp*3 + 0]), Z ero_one(pm2(L->Ph L do.)*L->d[i*3 + 


i])); 


r REMOVE "BIAS” REFERENCE H n*3+211V 

temp=min(temp, zero_one(L->x[2] L->d[i 3+2])), 


V 


return(temp); 


case 21 


temp 


lmin( zero_one(ps1 (L->Phi)*L->d[i*3 + 01), zero_one(ps2(L->Phi_dot)*L->d[i*3 + 


1 ])); 


r REMOVE "BIAS" REFERENCE ^ dn *3+211V 
temp=min(temp, zero_one(L->x[2] L->d[i 3+2])), 


V 


return(temp); 


“Imp-mint zero_one(ps1 (L->Phi)*L->dp*3 + 0]). zero_one(zo2<L->Phi_dot)*L->dP"3 + 


1])): 


r REMOVE "BIAS" REFERENCE H fi*34.2nv 

temp=min(temp, zero_one(L->x[2] L->dfi 3+2])), 


V 


return(temp); 


0386 23 min( zero_one(pm1(L->Phi)*L->dti*3 + 0]). zero_one( P b2(L->Phi_dotrL->dp*3 + 


temp* 


I]))*. 


r REMOVE "BIAS" REFERENCE 

temp-min(temp, zero_one(L->x[2] L->d[i 3+2])), 


7 


return(temp); 


0386 24: _one(pm1 (L->Phi)*L->dp*3+01), zero_one(pm2(L->Phi_dot)*L->d[r3 


temp=min( zero 

r REMOVE "BIAS" REFERENCE +21 ». 

temD=min(temp, zero_one(L->x[2] L->d[i 3+2])), 


retum(temp); 


0386 25: >_one(pm1 (L->Phi)*L->dp*3+0]), zero_one(ps2(L->Phi_dotrL->dP*3 + 


temp=min( zero. 


T REMOVE "BIAS" REFERENCE >drr3+2m - 
temD=min(temp, zero_one(L->x[2] L >d[ ])). 


7 


temp=min(temp 


return(temp); 


8386 261 >_one(pm1 (L->Phi)*L->dp*3 + 0]). zero_one(zo2(L->Phi_dot)*L->dP'3 + 


temp=min( zero. 


1])): 
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r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i 3+2])); 

7 

return(temp); 

C ^mp=min( zero_one(pb1 (L->Phi)*L->d[i*3+0]), zero_one(pb2(L->Phi_dot)*L->d[i*3+ 

1 ])): 

r REMOVE "BIAS” REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
case 28* 

temp=min( zero_one(pb1 (L->Phir L->d[r3+0]), zero_one(pm2(L->Phi_dot)*L->d[i*3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
case 29* 

temp=min( zero_one(pb1 (L->Phi)*L->dp*3+0]), zero_one(ps2(L->Phi_dot)*L->d[i # 3+ 

i])); 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 
caso 30* 

temp-min( zero_one(pb1 (L->Phi)*L->dp*3+01). zero_one(zo2(L->Phi_dot)*L->d[i*3+ 

1 ])): 

r REMOVE "BIAS" REFERENCE 

temp=min(temp, zero_one(L->x[2]*L->d[i*3+2])); 

7 

return(temp); 


double calculate_z_array(i, L) 
int i; 

LEARN_CYCLE * L; TIN: 7 
switch (i) { 

case 0i 

return(ps3(L->w[0])) ; 

case 1 ! 

return(ps3(L->w[1 ])) ; 
case 2 * 

retum(pm3(L->w[2])) ; 
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case 3: 

return(pm3(L->w[3])); 
CSS6 4 1 

return(ps3(L->w[4])) ; 
case 5i 

return(ps3(L->w[5])) ; 
css© 6: 

return(pm3(L->w[6])) ; 
case 7* 

return(pm3(L->w[7])) ; 
case 8 \ 

return(ps3(L->w[8])) ; 
case 9: 

return(ps3(L->w[9])); 
case 1 0: 

return(ps3(L->w[1 0])); 
case 111 

return(ps3(L->w[1 1])); 
case 1 2: 

return(nm3(L->w[1 2])) ; 

case 1 3: 

return(nm3(L->w[1 3])); 

pacp 1 4* 

return(ns3(L->w[1 4])) ; 
case 15: 

return(zo3(L->w[1 5])); 

poep 1 fi* 

retum(ps3(L->w[1 6])): 
case 1 7i 

return(pm3(L->w[1 7])) ; 
case 1 81 

return(pm3(L->w[1 8])); 
case 1 9*. 

return(ns3(L->w[1 9])); 


poep 20* 

return(ns3(L->w[20])): 

C rltum(ns3(L->w[21])); 

^IfurnlnsStL^wpg)); 
case 23: 

return(nm3(L->w[23])) ; 


case 24 1 

return(nm3(L->w[24])) ; 

case 25: 

return(ns3(L->w[25])); 


learn_cycl«.c 
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rocp PR - 

return(ns3(L->w[26])); 
C3S6 27 1 

return(nm3(L->w[27])) ; 
C3.S6 28* 

return(nm3(L->w[28])) ; 
case 29: 

return(ns3(L->w[29])); 
pqcp qry 

return(ns3(L->w[30])); 


} 

} 
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r IDENTIFICATION: Hana Shehadeh UnCon, Corporation*/ 

p June 1991 f 

r PURPOSE:*/ 

#define RAND_MAX 32767 

typedef struct { 


/•PARAMETER DECLARATION: 
/*TYPE VARIABLE[SIZE] 


/* 



int 

start.state; 

r 

int 

learn_flag; 

r 

int 

nl; 

r 

int 

n2; 

r 

double 

Phi; 

r deg 

7 

double 
* / 

Phi_dot; 

Tdeg 

/ 

double failure ; 

r 

double 

x_old[3] ; 

r 

double 

y new[31] ; 

r 

i* 

double 

zi[3l] ; 

r 

double 

ii [31] : 

r 

«* 

double 

el [31]; 

r 

Hi 

double 

w[31]; 

r 

fit 

double 

z[3l]; 

r 

fit 

double 

<[3i]; 

r 

double 

x[3]; 

r 

Hr - 

double a[93]; 

r < 

fit , 

double d[93]; 

r < 

fit 

double 

b[3]; 

r 

fit 

double 

c[3i]; 

r 

fit 

double 

e[3l]; 

r 

fit 

double 

y[3l]; 

r 

p 

double 

v_new ; 

/ 

double 

r_hat ; 

T 

double 

push ; 

r 

Ht 

double 

unusualness ; 

r 

p 

double 

sum ; 

/ 

fit 

double 

suml ; 

r 

fit 

double 

sum2 ; 

r 

fit 

double 

sum3 ; 

r 


****.....*.**.*~~*~*~*~****7 

USAGE BY DRIVER MODULES:*/ 
<INOUT> DESCRIPTION*/ 

7 


< > 

< > 

< > counter 

< > counter 

< > angle error 


7 

*/ 

7 

7 


< > rate error 



> 

> 

> 

> 

> 

> 

> 

> 

> 

> 


7 

7 

7 

7 

7 

*/ 

7 

7 

7 


7 

7 

> 

> 

> 

> 

> 

> FORCE 

> 

> 

> 

> 

> 


7 

7 

7 

7 

7 

7 

applied */ 

7 

7 

7 

7 


*/ 
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double sum4 ; 
double sumd ; 
double factorl ; 
double factor2 ; 
double factor3 ; 
double denom ; 
double v; 
double p; 
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r < > 

r < > 

r < > 

r < > 

r < > 

r < > 

r failure related parameter 7 
r failure related parameter 7 


} LEARN_CYCLE ; 



aat_md_wgta . c 


Mon Hay 4 10:23:00 1992 


1 


#include <stdio.h> 

#include <math.h> 

#include <sys/types.h> 

r EXTERNAL DATA STRUCTURE 

#include "../orb_fuzzy/learn_cycle.h" 


DEFINITION 7 


double md() /* Returns a floating-point between 0 and 1 , including 0. 7 

r rand is a number between 0 and 2 A 31 -17 
return ((double) rand() / (double)(RAND_MAX)); 

} 


set_rnd_weights (L) 

LEARN_CYCLE * L; 

{ 

int i,j; 

for (i = 0; i < 3; i++) 

L->b[i] = rnd() * 0.2 - 0.1 ; 
for (i = 0; i < 31 ; i++) 

L-xffl « md() * 0.2 - 0.1 ; 

L->e[i] = md() * 0.2 - 0.1 ; 

L->f[i] = md() * 0.5 + 0.5; 

L->w[i] = 1 .0; 

^ L->x_old[0]= L->x_old[1 ]= L->x_old[2]=0.0; 

for (j = 0; j <93; j++) 

{ 

L->aO] = rnd() * 0.2 - 0.1 ; 

L->dQ] = 1.0; 

} 


} 


Appendix B. Plots of Attitude-hold Test case 



SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


LVLH EULER PYR ROLL vs TIME 

RUN: Att Hold Vernier (light DB) 


-6.93889e-17 


LVLH -0.2 


















SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


LVLH EULER PYR PITCH vs TIME 

RUN: Alt Hold Vernier (Tight DB) 



VEHICLE: ORB_FUZZ_BATCH.*uie 
DATA SAMPLING FREQUENCY: 0100 Hi 


ORBITAL OPERATIONS SIMULATOR 


2 


Wed Apr 29 1992 04:35:26 PM 



SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


LVLH EULER PYR YAW vs TIME 

RUN: Alt Hold Vernier (Tight DB) 



TIME (minute*) 

VEHICLE: ORB_FUZZ_BATCHjuie 
DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


3 


Wed Apr 29 1992 04:35:26 PM 




SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


LVLH BIASED BODY ROLL RATE vs TIME 

RUN: Att Hold Vernier (Tight DB) 



TIME (minutes) 

VEHICLE: ORB_FUZZ_BATCH.ittte 
DATA SAMPLING FREQUENCY: 0.500 Hi 


ORBITAL OPERATIONS SIMULATOR 


4 


Wed Apr 29 1992 04:35:26 PM 



SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


LVLH BIASED BODY PITCH RATE vs TIME 

RUN: Alt Hold Vernier (Tight DB) 



TIME (minutes) 

VEHICLE: ORB_FUZZ_BATCRjutc 
DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


5 


Wed Apr 29 1992 04 35:26 PM 




SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


LVLH BIASED BODY YAW RATE vs TIME 

RUN: Au Hold Vernier (Tight DB) 



VEHICLE: ORB_FUZZ_BATCH-«uie 
DATA SAMPLING FREQUENCY: 0.500 Hr 


ORBITAL OPERATIONS SIMULATOR 


6 


Wed Apr 29 1992 04:35:26 PM 








SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


atter[2] vs TIME 

RUN: Att Hold Vernier (Tight DB) 



TIME (sec) 

MODULE: ORB_FUZZ_BATCHJ«w 
DATA SAMPLING FREQUENCY: 0.500 H* 


ORBITAL OPERATIONS SIMULATOR 


3 


Wed Apr 29 1992 04:35:44 PM 




SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


rter[0] vs TIME 

RUN: Alt Hold Vernier (Tight DB) 


2.03288c-19 
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0 100 200 300 400 500 600 700 800 900 1000 

TIME (tec) 


MODULE: ORB_FU2Z_BATCHJiw 
DATA SAMPLING FREQUENCY: 0.500 Hi 


ORBITAL OPERATIONS SIMULATOR 
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SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


rter[l] vs TIME 

RUN: Alt Hold Vernier (Tight DB) 



MODULE: ORBJFUE_BATCHiiw 
DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


5 


Wed Apr 29 1992 04:35:44 PM 




SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


rter[2] vs TIME 

RUN: Att Hold Vernier (Tight DB) 



MODULE: ORB_FUZZ_BATCHi«w 
DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


6 


Wed Apr 29 1992 04:35:44 PM 




00 


TIME (sec) 

MODULE: ORB_FUZZ_BATCHJe*m2 
DATA SAMPLING FREQUENCY: OJOO Hz 


ORBITAL OPERATIONS SIMULATOR 


h 


Wed Apr 29 1992 04:35:55 PM 










SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


w[15] vs TIME 

RUN: Att Hold Vernier (Tight DB) 



MODULE: ORB_FUZZ_BATCHJeam2 
DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


II 
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TIME (*ec) 

MODULE: ORB_FU2Z_BATCHJeam2 
DATA SAMPLING FREQUENCY: 0-500 Hz 


ORBITAL OPERATIONS SIMULATOR 
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SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


push vs TIME 

RUN: Att Hold Vernier (Tight DB) 



DATA SAMPLING FREQUENCY: 0.500 H* 


ORBITAL OPERATIONS SIMULATOR 


15 


Wed Apr 29 1992 04:35:55 PM 



SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


unusualness vs TIME 

RUN: Att Hold Vernier (Tight DB) 


-1.3877 8e-16 1 1 1 1 1 1 

- 0.1 



-03 



unusualness -0.5 — 



““ 





.1 

0 100 200 300 400 500 600 700 800 900 

TIME (sec) 

MODULE: ORB_FUZZ_B ATCHieam2 

DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


16 
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SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


failure vs TIME 

RUN: An Hold Vernier (Tight DB) 



TIME (tec) 

MODULE: ORB_FUZZ_BATCHietni2 
DATA SAMPLING FREQUENCY: 0300 Hz 


ORBITAL OPERATIONS SIMULATOR 


17 


Wed Apr 29 1992 04:35:55 PM 















SIMULATION APPLICATION: FUZZY ORBITER BATCH APPUCATION (9/14/90) 


d[26] vs TIME 

RUN: Att Hold Vernier (Tight DB) 



0 100 200 300 400 500 600 700 800 900 1000 

TIME (sec) 

MODULE: 0RB_FUZZ3ATCHJcam2 
DATA SAMPLING FREQUENCY: 0300 Hz 


ORBITAL OPERATIONS SIMULATOR 


23 


Wed Apr 29 1992 04:35:55 PM 










SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 

d[42] vs TIME 

RUN: Alt Hold Vernier (Tight DB) 

1.00415 — 



TIME (tec) 

MODULE: ORB_FUZZ_BATCHJe»m2 
DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


27 
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00 200 300 400 500 

TIME (fee) 


MODULE: ORB _FUZ£_BATCHle»m2 
DATA SAMPLING FREQUENCY: 0.500 Hx 


ORBITAL OPERATIONS SIMULATOR 
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SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


d[45] vs TIME 

RUN: Alt Hold Vernier (Tight DB) 
1.00415 



TIME (sec) 

MODULE: ORB_FUZZ_BATCHJeam2 
DATA SAMPLING FREQUENCY: 0.500 Hz 


ORBITAL OPERATIONS SIMULATOR 


30 
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SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


w[22] vs TIME 

RUN: Alt Hold Vernier (Tight DB) 



MODULE: ORB_FUZZ_BATCHie»m2 
DATA SAMPLING FREQUENCY: 0.300 Hz 


ORBITAL OPERATIONS SIMULATOR 


14 


Wed Apr 29 1992 04:35:55 PM 





SIMULATION APPLICATION: FUZZY ORBITER BATCH APPLICATION (9/14/90) 


f[14] vs TIME 

RUN: Alt Hold Vernier (Tight DB) 



MODULE: ORB_FUZZ_BATCHie»m2 
DATA SAMPLING FREQUENCY: 0.500 Hi 


ORBITAL OPERATIONS SIMULATOR 


3 


Wed Apr 29 1992 04:35:55 PM 





SIMULATION APPLICATION: FUZZY ORB ITER BATCH APPLICATION (9/14/90) 


f[15] VS TIME 

RUN: Att Hold Vender (Tight DB) 



TIME (sec) 

MODULE: ORB FUZ7. RATT H 
DATA SAMPLING FREQUENCY: 0 .500 Hz 


ORBITAL OPERATIONS SIMULATOR 


4 


Wed Apr 29 1992 04:35:55 PM 







